Detecting localized damage in cantilevered structures under nonstationary ambient excitations via Gabor spectral mode transmissibility functions

A method based on Gabor spectral mode transmissibility functions (GSMTFs) is proposed to detect local damage in a cantilevered structure under nonstationary ambient excitations. Gabor transformation and singular value decomposition are used to reduce the influences of other vibration modes on Gabor spectral mode transmissibility functions and process nonstationary structural responses, respectively. A new state characteristic based on the fundamental structure frequency is formulated on the basis of the GSMTFs, eventually leading to the development of a new damage indicator. The probability density functions of the damage indicator for healthy and damaged states can be estimated from the measured data, and the receiver operating characteristic (ROC) curve derived from these probability distributions and the corresponding area under the ROC curve (AUC) are used to determine the damage location. A six-degree-of-freedom system and a typical transmission tower are numerically studied, and the results show that the proposed method can estimate the structural damage location under nonstationary random loads. The proposed method is further validated with a planar frame in the laboratory, which exhibits multiple damage elements via random force hammer excitations. The results show that the AUC values computed for certain parts of the structure containing the damaged elements are greater than those for other parts of the structure, indicating the effectiveness of the proposed method. Moreover, the proposed method is compared with the dot product difference (DPD) index, and the results from the laboratory planar frame demonstrate that the proposed method can better identify damage.

www.nature.com/scientificreports/been successfully applied to truss structures and cantilevered beams 15,16 .By combining natural frequencies and mode shapes, several new objective functions have been proposed 17 , and optimization algorithms have been used to identify structural damage 18 .Worden et al. 19 proposed TFs as the damage metric from measured white noise excitations, and the saw-cut damage in a stringer was successfully identified.Simon Chesné et al. 20 employed TFs to detect and localize element stiffness loss in a four-degree-of-freedom (DoF) spring-mass system and a simulated cantilevered beam.TFs, among all the dynamic features reported, have shown promising sensitivity to structural damage.However, TF algorithms are subject to the effects of the excitation location, the selection of the analysis frequency bandwidth and the environmental noise effect.
Stationary white ambient excitation was assumed in the above studies.Wang et al. 21adopted power spectral density functions to evaluate TFs, and their results revealed that the identification error of the lower mode order is small under white noise excitations.Iván Gómez Araújo et al. 22 proposed an operational modal analysis method based on power spectral density transmissibility matrices.It was demonstrated that the TFs at the natural frequencies converge to a ratio of modal vibration amplitudes.Mao 23 adopted TFs at the fundamental frequency of a 3-DoF spring-mass system to identify structural damage under stationary excitations.TFs at the fundamental frequency of a structure have been shown to be sensitive to local damage.
Ambient excitations from vehicles and wind loads are nonstationary 24 .Jun Luo et al. 25 proposed the use of wavelet spectral transmissibility functions to address the damage detection problem of structures under these excitations with a damage indicator.Nevertheless, the study was confined to a shear structure, and the effect of measurement noise was not considered.Furthermore, the selection of wavelet basis functions for the wavelet transform is subjective.The Gabor transform has a high time-frequency resolution, with the product of the temporal resolution and frequency resolution equal to 0.5.It is a popular time-frequency analysis method, and selection criteria for shape parameters in the Gabor transform have been proposed 26 .Recently, the Gabor transform and corresponding Gabor spectral mode transmissibility functions (GSMTFs) have been successfully applied to mode shape identification under stationary white noise excitation, with singular value decomposition (SVD) used to reduce the influence of other vibration modes on the Gabor coefficient of modal responses at natural frequencies 26 .The error of identification for the lower mode order is smaller, and the GSMTFs at the natural frequencies converge to a ratio of amplitudes of vibration modes.
The Gabor transform and the GSMTFs at the fundamental frequency of the structure are introduced briefly for handling the nonstationary responses used to construct a new damage indicator for detecting damage in cantilever structures.The Gabor transformation and SVD are used to process the nonstationary dynamic responses and reduce the influence of other vibration modes on the Gabor coefficient of responses at natural frequencies 26 .The existing TF calculation methods are mostly based on power spectrum estimation, which cannot be used for nonstationary responses in theory.As a commonly used nonstationary signal processing method, the Gabor transform can avoid these disadvantages.The damage indicator is derived on the basis of the GSMTFs at the fundamental frequency.The ROC curve and the area under the ROC curve (AUC) are used to localize the damage.Two numerical examples and a laboratory structure are assessed to demonstrate the feasibility of the proposed algorithm.
The rest of this paper is organized as follows.Section "Basic theory of the damage detection method based on GSMTFs" presents the basic theory of the damage detection method, which is based on GSMTFs, followed by the derivation of the new damage indicator.Section "Numerical validation" validates the proposed approach with two numerical examples.Section "Laboratory verification" presents the experimental validation of the proposed approach with a laboratory frame structure.The results are compared with those obtained from the dot product difference (DPD) index.Conclusions are presented in Section "Conclusions".

Purification of GSMTFs via the Gabor transform and SVD
For a given signal y(t), its Gabor coefficients can be expressed as 27 where * denotes the complex conjugate, g is the basis function, b is the time shift, a is the frequency shift, N is the length of the window function and M is the number of discrete frequency points.A Gaussian window function is usually used as the basis function.The time-domain and frequency-domain expressions of the Gaussian window function can be expressed as where σ is the shape parameter of the Gaussian window function, which affects the frequency resolution and temporal resolution of the basis function and can be determined from the literature 26 .
For an n-DoF structure, the Gabor coefficients of the acceleration response at the i-th DOF, i.e., signal ÿ i , can be expressed as where ω r is the r-th modal frequency, ϕ ir is the modal amplitude of the r-th mode, m r and ξ r are the modal mass and modal damping ratio of the r-th mode, a r is the frequency shift corresponding to the r-th modal frequency, and ω dr is the r-th damped modal frequency.The r-th modal response is usually assumed to be dominated by that from the corresponding vibration mode with the approaches in Eq. ( 3).Therefore, the Gabor coefficients of signal ÿ i at the r-th modal frequency can be expressed as where G ÿi,r (b, a r ) and �G ÿi,r (b, a r ) are the Gabor coefficients of the r-th vibration mode and other vibration modes, respectively.
Furthermore, to reduce the influence of other vibration modes in the r-th modal response, a Gabor coefficient matrix H(a r ) is constructed, and SVD 28,29 is used to extract the purified matrix Ĥ(a r ) 26 , which is related only to the r-th modal response using the first left and right singular vectors, as shown below.
where Ĥ and ΔĜ denote the singular value decomposed Gabor coefficient matrices of the r-th vibration mode and other vibration modes of signal ÿ i , respectively.s i , u i and v i are the i-th singular values, i-th left singular column vectors, and i-th right singular column vectors of matrix H(a r ), respectively.L is the length of signal ÿ i .b e (e = 1,2,…,L) denotes the e-th time shift.Hence, the GSMTFs between the i-th and j-th DoFs with a k-th reference DoF can be defined via the purified Gabor coefficients of the modal responses in matrix Ĥ(a r ) as 26 where * denotes the conjugate.GSMTF k i,i-1 can be calculated via the purified Gabor coefficients of the modal responses ÿ i,r , ÿ j,r and ÿ k,r extracted from the purified Gabor coefficient matrices in Eq. ( 5).The purified Gabor coefficient matrices are obtained from all measured responses of the whole structure.As shown in Eq. (6), the GSMTF k i,i-1 (a r ) is equal to mode shape coefficient ratios ϕ ir /ϕ jr in theory.However, if the lower and higher mode contributions are considered, because of the influence of �G ÿi,r (b, a r ) , the GSMTF k i,i-1 (a r ) cannot be derived into the ratios ϕ ir /ϕ jr in theory, and the established damage indicator does not have a clear physical meaning.

The proposed damage indicator
After the structural natural frequencies are identified, the purified GSMTFs at the fundamental frequency are selected to form the damage indicator because they can be identified more accurately under ambient excitations.The reason for using the GSMTFs at the fundamental frequency is that there is a corresponding relationship between the proposed D i,i+1 (i = 1,2,…,n-1) and the curvature value of the mode shape, and the calculation of (3) www.nature.com/scientificreports/curvature values is prone to significant errors at the location where the mode coefficient symbol changes, which is more likely to appear in the second and upper modal orders.Taking a structure with n elements, as shown in Fig. 1, the state characteristics associated with the modal curvature of two adjacent connected elements can be defined as where k is the reference DoF.Any measured DOF in a structure can be used as the reference DoF, and the mean of GSMTF k i+1,i-1 of the Gabor spectral mode transmissibility functions between two DoFs can be calculated; subsequently, the state characteristics are D i,i+1 .As shown in Eq. ( 7), there is a corresponding relationship between D i,i+1 (i = 1,2,…,n-1) and the first-order mode shape at nodes in elements i and i + 1, as shown below.D i,i+1 represents the curvature value of the first-order mode shape at nodes in elements i and i + 1.Therefore, D i,i+1 can be used to characterize the healthy state of elements i and i + 1, and the changes in the state characteristics D i,i+1 under the initial state and damage state can be used to identify damage.
Since there are errors from measurement and computation, the calculated state characteristics D i,i+1 vary, even if the structure is not damaged.To reduce the influence of the above errors on damage identification, three structural states are defined, i.e., the initial state, healthy state and damaged state.The corresponding state characteristics can be expressed as D r i,i+1 , D u i,i+1 and D d i,i+1 .D r i,i+1 and D u i,i+1 are both calculated from the responses under undamaged structures.If Q samples of responses of the structure under initial, healthy and damaged states are collected, Q samples of indices D r i,i+1 , D u i,i+1 and D d i,i+1 will be obtained.Then, the mean value of D r i,i+1 can be calculated and used as the final identified index under the initial state, which can be expressed as Ď r i,i+1.Finally, the damage indicators C u i,i+1 and C d i,i+1 of the Q samples can be calculated as follows: where C u i,i+1 and C d i,i+1 are the damage indicators for the pair of elements i and i + 1 under healthy and damaged states, respectively.GSMTF, D and C are matrices or vectors, and GSMTF k i,i-1 , D i,i+1 and C i,i+1 are scalers.Furthermore, the probability density functions (PDFs) of the damage indicators C u i,i+1 and C d i,i+1 are assumed to be normally distributed as  Therefore, the proposed damage detection method is feasible because it checks whether the PDFs of C u i,i+1 and C d i,i+1 have been changed, as noted from the ROC curve and AUC values.The ROC curve 30 can then be drawn from the PDFs for the healthy and damaged states of the structure, and the AUC is used to determine the damage location.The AUC has a maximum value of unity, which corresponds to a complete damage state, and a value of 0.5 corresponds to an undamaged state of the pair of elements considered.More detailed discussions on the theory of AUCs and ROCs can be found in 30 .
The flowchart of damage diagnosis of the proposed algorithm is shown in Fig. 1 and described below.
1.The index D r i,i+1 for the intact structure is calculated from the purified GSMTFs.Q samples of the responses of the structure under nonstationary excitations are collected, and the recommended value of Q should be greater than 33 for good curve fitting 31 of the damage indicator PDFs in Eq. ( 9).The value of Q is determined on the basis of the total length of the data and should be as large as possible.For each sample of responses, the Gabor coefficients of the signals are calculated with the matrix H(a r ) obtained from Eq. ( 5).The purified GSMTFs and the state characteristics D r i,i+1 are obtained from Eqs. ( 6) and ( 7

Numerical validation A 6-DoF spring-mass system
Damage scenarios studied and nonstationary excitations A 6-DoF spring-mass system, as shown in Fig. 2, is studied.The weights of all masses are equal to 1.0 kg, and the spring stiffnesses are all equal to 1500 N/m.Proportional damping is assumed with , α = 0.3081 and β = 7.5 × 10 -4 .Local damage to the structure is simulated by reducing the spring stiffness values, and the five damage scenarios studied are shown in Fig. 2.
A horizontal nonstationary force is applied to mass 6 without loss of generality.It is generated from the filtering of white noise via an amplitude modulation function as follows: The horizontal acceleration responses at all the DoFs under the applied force are calculated via the Newmark method with a sampling frequency of 50 Hz for a period of 1000 s.Fifty sets of structural responses for the intact and healthy states of the structure and when the structure has one of the damage scenarios shown in Fig. 2 are similarly obtained.
The response of mass 6 without noise is shown in Fig. 3 together with its power spectral diagram for illustration.The fundamental frequency of the system is noted at 1.466 Hz.The power spectrum is estimated via the modified periodogram method.The number of discrete Fourier transform points is 1024, and the Hanning window is used.This signal processing technique is adopted in all the studies in this paper.

Damage detection without noise
The number of discrete points in the Gabor transform is 1024, with a shape parameter of 2.306.The purified GSMTFs at the fundamental frequency without noise are calculated from the intact, healthy and damage states  6) and ( 7).The PDFs of C u i,i+1 and C d i,i+1 for damage case 1 are calculated from Eqs. ( 8) and ( 9), as shown in Fig. 4. The distributions of the two indicators are significantly different for partial structures 1 and 2. This is because spring k 2 within partial structures 1 and 2 is damaged in this scenario.
The ROC curve of the PDFs of the damage indicators for damage case 1 is shown in Fig. 5, with the vertical axis representing the true positive rate and the horizontal axis representing the false-positive rate 30,32 .It is noted 33   The ROC curves of partial structures 1 and 2 in Fig. 5 almost overlap with the upper left corner, and those for other partial structures vary around the diagonal line.The AUC values from these curves are calculated and shown in Table 1 and Fig. 6.The AUC values of the partial structures, including the damage elements, are much greater than 0.85 32 , indicating the presence of damage in these partial structures.Similar observations are noted in the results of all damage scenarios studied, with AUC values associated with damage approaching unity, whereas other values not related to any damage are approximately 0.   11).The ratio of the noise amplitude to the simulated acceleration (NSA) is set to 0.10, 0.20 and 0.30 for this study 26 .
where δ is the NSA and max(•) is the absolute maximum value of the signal.randn is a sample from a normal random distribution.
The proposed method is conducted with noise-polluted responses, and the PDFs of damage indicators C u 1,2 and C d 1,2 for damage scenario 1 are shown in Fig. 7 for illustration.The AUC values for different damage scenarios are shown in Fig. 8.The red horizontal plane corresponds to the threshold value of 0.85 for damage occurrence.The values for the partial structures that include the damage elements are much larger than those without any damage, and the AUC values of the damaged partial structures decrease with increasing noise level.For example, the stiffness k 5 is damaged in scenario 3, which affects partial structures 4 and 5.When the noise level increases from 0.1 to 0.3, the AUC values of partial structure 4 decrease from 0.996 to 0.787.The AUC value is above the threshold value of 0.85 when NSA < 0.2, indicating the robustness of the proposed method to measurement noise interference.

An electricity transmission tower
The structure and damage scenarios studied A 2B3-ZM2 electricity transmission tower, as shown in Fig. 9, is taken as a typical large-scale structure to verify the effectiveness of the proposed method.The structure is 5 × 5 m in length at the bottom, with a total height of 21 m.It is a lattice structure consisting of many vertical, horizontal and inclined members.Sensors are placed at 8 levels horizontally to collect the dynamic responses under horizontal loading.The structure is divided into 7 partial structures, as shown.Notably, in this example, the proposed method is not applicable for detecting local damage above sensor S8.Local damage to the structure is simulated by reducing the elastic modulus of certain (11)

Nonstationary excitation and dynamic responses of the tower
To obtain the responses of the tower under nonstationary excitation, a nonstationary fluctuating wind load is assumed, which is obtained via the Kaimal wind spectrum and a normal turbulent wind model based on the IEC61400-1 (Edition 3) code 34 .The turbulence characteristics are set as Class A to simulate a larger fluctuation.
The assumed wind load is applied along the x-direction to the two top joints of the tower at the level of sensor S8.The acceleration responses from all 8 sensors along the x-direction are obtained from a finite element  www.nature.com/scientificreports/model of the structure.All the structural members between the ground level and sensor 2 are included in partial structure 1.For the i-th partial structure, all the horizontal and inclined members between sensors i-1 and i + 1 are included.One hundred sets of data from the healthy state and 50 sets from each damage state are studied.These sets of responses are used to extract the purified GSMTFs and index D i,i+1 under different states.Furthermore, to consider the influence of measurement noise, three cases with polluted measurements from different noise levels of 0.10, 0.20 and 0.30 are also studied.The sampling frequency is 50 Hz, and the sampling duration is 100 s.The assumed wind load acting at the level of sensor S8 is shown in Fig. 10(a), and the acceleration response at S8 in the healthy state is shown in Fig. 10(b).The response is noted to exhibit amplitude nonstationarity, with notably larger amplitudes in the middle of the time period.
The power spectral diagram of the acceleration response is shown in Fig. 10(c), which shows the fundamental frequency of the system at 2.444 Hz.

Damage detection results
The Gabor transform on the acceleration responses is conducted similar to those in previous studies.Fifty samples of structural responses without noise are simulated for the initial, healthy or damage states of the structure.
With the proposed method, the purified GSMTFs at the fundamental frequency of the structure under initial, healthy and damage states can be determined, and 50 sets of state characteristics D r i,i+1 , D u i,i+1 , and D d i,i+1 can be calculated.The ROC curve and the AUC values are calculated similarly to those shown in the last example.
The calculated AUC values are shown in Table 2 and Fig. 11.The partial structure containing the damage elements has much greater AUC values than those without the damage elements.Especially for cases without noise, the AUC values of partial structures 3 and 4 are very close to the maximum value of unity for all damage cases, indicating that the elements in these partial structures are damaged.The identified results are consistent with the simulated damage location.With a 10% noise level, only the damaged partial structures under damage case 2 can be identified, and the lower degree damage under damage case 1 cannot be identified.However, with The results show that noise has an impact on the damage diagnosis of large actual structures.The reason is that the real transmission tower is a lattice structure, and the proposed method in this paper assumes that the cantilevered structure consists of 8 partial structures.For each partial structure, all the vertical, horizontal and inclined members between two adjacent sensors are included, and only the damage to the main chords of the transmission tower is simulated.Therefore, the degree of decrease in the actual lateral stiffness is less than the preset damage degree, and the identification results are more susceptible to noise influence than the 6-DoF spring-mass system is.

Steel frame testbed and damage cases
A planar laboratory steel frame, as shown in Fig. 12, is studied to verify the proposed method.The frame is fabricated from steel plates, gusset plates and bolts.The dimensions of all the vertical steel members are 350 × 65 × 4 mm.
Three damage cases are simulated in the test by replacing the intact vertical member with one whose width is reduced along its full length.Damage member 1 has the width of the member reduced by 20% on one edge, whereas damage member 2 has the width reduced by 20% on both edges.The structure is divided into three partial structures, as shown for damage identification.

Excitation and dynamic responses
Random taps by two hammers at the first and second floor levels of the frame are used to excite the structure at approximately 0.5 s intervals.The material at the hammer head is soft rubber.The acceleration responses of each floor were measured by YD81D-V accelerometers installed at the joints and collected by a dynamic data acquisition system.The sampling frequency is 150 Hz, and the total sampling duration is 3200 s for measurement in the healthy state.These responses are divided into 80 segments of 40 s each for the estimation of the PDFs of the damage indicators.
The acceleration responses of sensors S1 and S4 are shown in Fig. 13(a),(b).The power spectrum of the acceleration response at sensor S4 is shown in Fig. 13(c), which shows the fundamental frequency of the system at 4.470 Hz.It is similar to those in the previous study but has 2048 discrete points in the spectrum from the 6000 data points in the data segment of the signal.The first 40 segments are used to calculate the purified GSMTFs at the fundamental frequency of the structure and the state characteristics D r i,i+1 in the initial state.The remaining 40 segments are used to calculate the purified GSMTFs at the fundamental frequency and state characteristics D u i,i+1 in the healthy state.Then, the PDF of damage indicators C u i,i+1 in the healthy state can be identified from Eqs. ( 8) and ( 9).There are 2048 discrete points in the Gabor transform, with the shape parameter of the Gabor transform equal to 1.537.
For each damage scenario studied, the total sampling time is 1520 s, and 38 data segments are obtained with 40 s each.All 38 segments are used to calculate the purified GSMTFs at the fundamental frequency and state characteristics D d i,i+1 for each damage scenario.Then, the PDF of damage indicators C d i,i+1 for each damage case can be estimated on the basis of the state characteristics D r i,i+1 and D d i,i+1 .The ROC curves are then plotted from the PDFs of damage indicators C u i,i+1 and C d i,i+1 , and the AUC values calculated for each damage scenario are shown in Table 3 and Fig. 14.The results show that the AUC values of the partial structure containing the damage elements are larger than those without the damage elements.In damage cases 1 and 2, only the AUC values of the partial structures containing the damage elements are greater than 0.85, which is consistent with the true location of the damage.For damage case 3, only the AUC values of the partial structure containing the damage element with 40% width reduction are greater than 0.85, whereas those with 20% width reduction at the third floor cannot be identified.

Comparison with existing dot product difference indices
The dot product difference (DPD) index has also been proposed to identify structural damage on the basis of GSMTFs or TFs at the fundamental frequency 23 .The index has the advantage of scaling between 0 and 1.The DPD index for damage identification can be expressed as Eq. ( 12) when GSMTFs at the fundamental frequency are used.
where T i,j is the GSMTF at the fundamental frequency between the i-th and j-th DoFs, which is equal to the mean value of GSMTF k i,j (a 1 ) (k = 1,2,3,4).Superscripts r and d denote the initial state and damage state, respectively.Similarly, the GSMTFs at the fundamental frequency under initial and healthy states can be used to reduce the influence of noise.The DPD index under the healthy state can be expressed as Eq. ( 13).( 12)  4 and Fig. 15.The results show that the AUC values for all partial structures are less than 0.85, making the identification of damage locations difficult.Therefore, the proposed method can better identify damage.

Conclusions
In this paper, a new damage detection method based on the Gabor transform, SVD and ROC curves is developed.A new state characteristic is established on the basis of the purified GSMTFs at the first natural frequency, and a damage indicator is proposed.Furthermore, the PDFs of damage indicators under healthy and damaged states are identified, and the ROC curves and AUC values are used to determine the damage location.Two numerical examples and a plane frame in the laboratory are used to verify the effectiveness of the proposed method.
1.In the six-degree-of-freedom system numerical example, amplitude nonstationary loads are applied.Different damage locations and degrees are considered, and the influence of measurement noise is studied.The results show that the proposed method can identify the damage location and is robust to noise.Even at the 30% noise level, the AUC values of the damaged partial structures are still much greater than those of the other structures.2. In the transmission tower numerical example, the assumed nonstationary fluctuating wind loads and two damage cases are considered.The results show that the damage locations can be identified successfully.Moreover, because of the simplification of the identification structure, damage identification for minor damage to local components is more susceptible to noise interference, so measurement equipment with an antinoise system should be adopted.3.In the laboratory test, random percussion excitation is considered, and the damage location and degree are studied.The results show that the AUC values of the partial structure containing the damage elements are much greater than those of the other structures are, and the proposed method can identify the damage location.Moreover, the proposed method is compared with the DPD index, and the results show that the AUC values for all partial structures are less than 0.85, making the identification of damage locations difficult.The proposed method can better identify damage.www.nature.com/scientificreports/However, in the proposed method, the structures are simplified as cantilever structures with several DoFs, and only the damage between three adjacent DoFs can be considered, not for each member.Therefore, the applicability of the proposed method for identifying other types of structures or damage members requires further in-depth research.

Figure 1 .
Figure 1.Flowchart of damage diagnosis using purified GSMTFs of the fundamental frequency.

Figure 3 .
Figure 3. Acceleration response and power spectrum of mass 6 without noise.

Figure 4 .
Figure 4. PDFs of the damage indicators under damage case 1.

Figure 5 .
Figure 5. ROC curves of the probability distributions of the damage indicators for damage case 1.

Figure 7 .
Figure 7.The PDFs of indicators C u 1,2 and C d 1,2 for damage scenario 1.

Figure 8 .
Figure 8. AUC values for different damage scenarios with noise.

Figure 10 .
Figure 10.One of the simulated wind loads and acceleration responses S8 under the healthy state.

Figure 11 .
Figure 11.AUC values for different damage cases with different noise levels.

Figure 12 .
Figure 12.The laboratory steel frame and damage cases.

( a )Figure 13 .
Figure 13.Acceleration responses and their power spectra in the healthy state.

,Table 4 .Figure 15 .
Figure 15.AUC values for different damage cases according to the DPD index.
where subscripts u and d denote the healthy state and damaged state, respectively.μ u and p u are the mean and standard deviation, respectively, of C u i,i+1 , where μ d and p d are the mean and standard deviation of C d i,i+1 , respectively.

Table 1 .
AUC values for different damage cases.Significant values are in bold.

Table 2 .
AUC values for different damage cases with different noise levels.Significant values are in bold.

Structure Floor AUC values Damage Case 1 Damage Case 2 Damage Case 3
the superscript u denotes the healthy state.Similarly, if Q samples of responses of the structure under initial, healthy and damaged states are collected, Q samples T r i,j , T u i,j , and T d i,j can be calculated.Then, Q samples of DPD u i,j and DPD d i,j can be obtained.Finally, the ROC curve and AUC values can be calculated.The mean values of GSMTF k 2,1 , GSMTF k 3,1 , and GSMTF k 4,2 (k = 1,2,3,4) under different structural states are used to calculate the DPD index and AUC values, as shown in Table 3. AUC values for different damage cases.Significant values are in bold.Figure 14.AUC values for different damage cases.where